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Abstract 

Background: The protein secretory pathway must process a wide assortment of native proteins for eukaryotic cells 
to function. As well, recombinant protein secretion is used extensively to produce many biologies and industrial 
enzymes. Therefore, secretory pathway dysfunction can be highly detrimental to the cell and can drastically inhibit 
product titers in biochemical production. Because the secretory pathway is a highly-integrated, multi-organelle 
system, dysfunction can happen at many levels and dissecting the root cause can be challenging. In this study, we 
apply a systems biology approach to analyze secretory pathway dysfunctions resulting from heterologous 
production of a small protein (insulin precursor) or a larger protein (a-amylase). 

Results: /-//\C7-dependent and independent dysfunctions and cellular responses were apparent across multiple 
datasets. In particular, processes involving (a) degradation of protein/recycling amino acids, (b) overall transcription/ 
translation repression, and (c) oxidative stress were broadly associated with secretory stress. 

Conclusions: Apparent runaway oxidative stress due to radical production observed here and elsewhere can be 
explained by a futile cycle of disulfide formation and breaking that consumes reduced glutathione and produces 
reactive oxygen species. The futile cycle is dominating when protein folding rates are low relative to disulfide bond 
formation rates. While not strictly conclusive with the present data, this insight does provide a molecular 
interpretation to an, until now, largely empirical understanding of optimizing heterologous protein secretion. This 
molecular insight has direct implications on engineering a broad range of recombinant proteins for secretion and 
provides potential hypotheses for the root causes of several secretory-associated diseases. 
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Background 

The protein secretory pathway is an extensive process in 
eukaryal cells, as it is responsible for processing approxi- 
mately one-third of all proteins. Substantial cellular 
resources are therefore utilized to maintain this pathways 
functions, and stressed conditions in the secretory path- 
way have consequences for the whole cell [1]. Distress in 
secretory pathway organelles has been implicated as the 
molecular basis for several diseases, for example, (3 cell 
apoptosis in diabetes, cystic fibrosis, and prion-related 
disease, among others [2], In biotechnology, efficient 
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secretion of useful recombinant proteins in yeast and 
fungi is a key industrial objective with applications in 
enzyme production required for the production of bio- 
fuels, detergents, fabrics, food, and biologies, such as imu- 
noglobulins, hormones, and vaccines. Significant effort has 
gone into engineering yeast for increasing protein secre- 
tion [3]. Strategies, such as changing environmental para- 
meters (for example, temperature, media composition) [4] 
or altering genetics, can increase secretion for some pro- 
teins, but they rarely represent generic solutions for 
improving protein secretion [5,6]. The lack of a single 
engineering strategy that improves protein secretion across 
the board implies that there are several possible bottle- 
necks in the secretory pathway, and different proteins may 
be constrained in different ways. There is therefore a 
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requirement for more fundamental insight into this com- 
plex pathway that involves a very large number of 
components. 

In yeast, the secretory pathway is a multi-organelle sys- 
tem that is responsible for trafficking proteins to the 
extracellular space, cell membrane, or vacuole [7]. During 
this transit, multiple processes must be coordinated, 
including folding, specific proteolytic cleavage, glycosyla- 
tion, and disulfide bond formation, all with a layer of 
quality control at key check points. The pathway requires 
substantial cellular resources to perform these tasks, such 
as glycans, electron acceptors, electron donors, and ATP. 
In the ER, the nascent peptide is folded into its native 
structure while disulfide bonds are formed. The rate of 
protein folding is dependent upon the complexity of the 
protein to be folded, the availability of chaperones to 
assist folding, and ATP used by the chaperones [1]. Pro- 
teins that are slow to fold or terminally misfolded pro- 
teins are removed from the ER via the ER-associated 
degradation (ERAD) pathway [8]. Disulfide bond forma- 
tion requires the removal of electrons from cysteine 
thiols via protein disulfide isomerase (PDI) and Erolp to 
the final electron acceptor, typically oxygen [9,10]. This 
process produces reactive oxygen species (ROS) in stoi- 
chiometric amounts to the number of disulfide bonds 
formed [11]. Disulfide bond formation is random, and 
incorrect bond pairs must be exchanged for native bonds 
via PDI-based processes [12]. In addition, reduced glu- 
tathione (GSH) acts as a buffer for the redox state of the 
ER [13]. A more detailed description of oxidative protein 
folding can be found in the reviews by Sevier et al. and 
Chakravarthi et al. [14,15]. 

The secretory pathway must adjust the chaperone 
capacity, oxidizing equivalents, ATP, glycan, and other 
metabolic requirements, as well as trafficking patterns, 
based on the portfolio of proteins that need to be 
expressed at a given time, and the resources required to 
process that set of proteins. In yeast, the unfolded protein 
response (UPR) is one transcriptional mechanism that 
adjusts secretory resources and controls to handle over- 
load of the folding machinery in the ER [16]. In the UPR, 
accumulation of unfolded proteins in the ER signals a 
pathway that results in translation of Haclp, a transcrip- 
tion factor (TF) known to activate or repress over 100 
genes, including many ER-associated proteins such as 
Kar2p, Pdilp, and Erolp [17]. 

In this study, we identified biological mechanisms 
which alter the secretory pathway in response to secre- 
tion of recombinant proteins with different properties 
(size, number of disulfide bonds, and glycans) in a 
Haclp-dependent and independent manner. The secre- 
tory pathway was perturbed by secreting a small protein, 
human insulin precursor (IP), or a comparatively larger 
protein, a-amylase, in wild-type (WT) and Ahacl 



Saccharomyces cerevisiae. These proteins were chosen 
because the two proteins elicit different behavior in the 
secretory pathway. These differences will arise because 
a-amylase is a relatively larger (and likely more difficult 
to fold), has an odd number of cysteines (which may 
complicate disulfide isomerization) and has glycosylation, 
compared to insulin which is small, has even number of 
cysteines, and is not glycosylated. As well, a-amylase has 
one more disulfide bond than IP. To identify biological 
mechanisms, we characterized changes in physiological 
properties (specific growth rate, carbon utilization effi- 
ciency, and recombinant protein secretion), TF activity 
(as inferred from transcriptome analysis) and metabolic 
demand (as inferred by changes in metabolic flux diver- 
sion). Through this, we identified the following biological 
processes: amino acids recycling from degraded proteins, 
trans-Golgi network (TGN) sorting changes, overall 
expression repression, and oxidative stress. Motivated by 
secretory-related oxidative stress observations, we pre- 
sent a model for disulfide bond formation and electron 
transfer in the ER which takes into account thermody- 
namic irreversibilities caused by differences in electron 
affinity. The proposed model explains the non-stoichio- 
metric ROS formation that we observed that results from 
disulfide bond formation and causes oxidative stress 
under folding-stress conditions. If proven by genetic and 
biochemical results, the futile cycle model yields insight 
into a fundamental problem in secretory stress and 
reveals new avenues to reduce oxidative stress and 
increase productivity in industrial protein production. 

Results 

Protein size and Had p activity affect protein secretion 
quantity and cell growth 

Yeast strains were constructed that produce and secrete 
(a) IP or (b) a-amylase and were compared to yeast 
strains containing (c) an empty vector in both wild-type 
and HAC1 deletion backgrounds. IP and a-amylase were 
chosen because they are very different types of proteins 
to secrete. IP is 51 amino acids in length, with six 
cysteines forming three disulfide bonds, and no glycosyla- 
tion. a-amylase is 478 amino acid in length, with nine 
cysteines forming only four disulfide bonds and one gly- 
cosylation. The odd number of cysteines in a-amylase 
complicates disulfide pairing, as the random isomeriza- 
tion process may incorporate the cysteine that should 
not be incorporated into a disulfide bond. Both proteins 
were targeted for secretion using a YAP3 pre sequence 
(21 amino acids, cleaved off in the ER) and a rationally 
designed pro sequence (TA57, 42 amino acids, no glyco- 
sylation or disulfides) were cloned behind a TDH3 pro- 
moter in a high copy 2 micron plasmid [18]. a-amylase 
was expressed using the same plasmid, promoter, and 
leader sequences. These strains are named WN (WT 
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with empty vector), Wl (WT secreting IP), WA (WT 
secreting a-amylase), dN (Ahacl with empty vector), dl 
(Ahacl secreting IP), and dA (Ahacl secreting a-amy- 
lase). Strains were characterized in batch fermentation to 
understand the effects on cell physiology. 

The cellular burden induced by (a) synthesizing and 
secreting IP and a-amylase and (b) deleting the key TF 
for the UPR, Haclp, substantially affected the cells. Pro- 
tein titers in WT strain were 9 mg/L and 20 mg/L, for IP 
and a-amylase, respectively (Figure la). On a per biomass 
basis, this is approximately half the insulin produced, and 
one-third the a-amylase reported for rich media [19,20]. 
Rich media appears to be favorable for heterologous pro- 
tein production, but may present complications in down- 
stream separations. Comparing the small and larger 
proteins, a-amylase was secreted in higher levels on a 
mass basis, but six-fold more insulin molecules were 
secreted (1.52 uM IP in Wl compared to 0.26 uM a-amy- 
lase in WA). Ahacl strains secreted significantly less pro- 
tein than WT, confirming that Haclp is important for 
efficient secretion (Figure la) [5]. 

Reduced specific growth rates imply impairment of cel- 
lular processes (Figure lb). In WT yeast, IP production 
did not affect growth; however, a-amylase production 
reduced growth by 25%. This, combined with the differ- 
ences in protein titers, implies that a-amylase is more 
challenging to fold and secrete than IP. In the Ahacl 
background, recombinant protein strains dl and dA had 
approximately 20% lower growth rates compared to dN. 
This growth reduction occurs despite no change in speci- 
fic glucose uptake rate (Additional file 1, Tables SI and 
S2) pointing toward higher energy requirements to main- 
tain homeostasis in Ahacl while trying to secrete recom- 
binant proteins. Ahacl strains had overall lower final cell 



densities. Ahacl strains produced more glycerol than 
WT strains implying impaired oxidative processes in the 
Ahacl strains (Additional file 2). 

Secretory stress shifts metabolism to increase oxygen and 
ATP requirements 

The physiological changes due to the secretory pertur- 
bations affect the distribution of resources through the 
metabolic network. The glucose uptake and range of 
products produced were altered by the protein produc- 
tion conditions (Table 1). Changes in the underlying 
metabolic network were estimated by flux balance analy- 
sis (FBA) using a yeast central carbon metabolism 
model, constrained by measured extracellular fluxes 
(Additional file 1, Tables SI and S2, Additional files 3 
and 4) [21]. Figure 2a shows a metabolic map of central 
carbon metabolism for each of the six conditions based 
on the exchange fluxes in Table 1 and the FBA analysis. 
The shift in metabolic fluxes were correlated with 
changes in redox requirements. As expected, the cata- 
bolic functions of the TCA cycle was predicted to have 
very low activity due to glucose repression [22]. Figure 
2b shows that the oxygen uptake was twice as high in 
the strains that were growth inhibited (for example, 
WA, dl, dA) than those that were not. This increased 
oxygen uptake was not used for oxidative phosphoryla- 
tion, as the biomass yields on glucose were lower in 
WA, dl, and dA, and it may therefore be a result of 
increased oxidation in connection with formation of dis- 
ulfide bonds. 

Figure 2c shows that the maintenance ATP consump- 
tion is increased in WA, dl, and dA according to FBA 
calculations. In WT background, Wl did not consume a 
detectable increase in ATP, likely because IP is short 




WN Wl WA dN dl dA 



WN Wl WA dN dl dA 
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Figure 1 Secretory perturbations affect yeast physiology (a) Final recombinant protein titer. Ahad strains were severely inhibited in 
recombinant secretion, (b) Specific growth rate on glucose. The combination of Ahacl and recombinant secretion had the most severe effect 
on growth, however even in wild-type background, a-amylase hindered growth, (c) Final cell concentration. Wild-type, no protein secretion 
(WN), wild-type insulin precursor secretion (Wl), wild-type a-amylase secretion (WA), Ahacl no protein secretion (dN), Ahacl insulin precursor 
secretion (dl), Ahacl a-amylase secretion (dA). Measurements are mean +/- s.e.m. (n = 3). 
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Table 1 Physiological parameters of recombinant protein secretion strains 3 



Strains 13 


M max [h" 1 ] 


Y S x 




Y S E 




Y SG 




Y S A 




Ysco2 




Carbon balance 


WN 


0.43 +/- 


0.014 


0.14 +/- 


0.001 


0.32 +/- 


0.041 


0.067 +/- 


0.009 


0.048 +/- 


0.0005 


0.30 +/- 


0.019 


0.89 


Wl 


0.40 +/- 


0.012 


0.13 +/- 


0.002 


0.35 +/- 


0.029 


0.055 +/- 


0.005 


0.056 +/- 


0.0047 


0.30 +/- 


0.014 


0.92 


WA 


0.32 +/- 


0.007 


0.11 +/- 


0.003 


0.31 +/- 


0.006 


0.060 +/- 


0.007 


0.049 +/- 


0.0023 


0.30 +/- 


0.002 


0.84 


dN 


0.38 +/- 


0.005 


0.13 +/- 


0.004 


0.37 +/- 


0.025 


0.046 +/- 


0.003 


0.035 +/- 


0.0046 


0.29 +/- 


0.020 


0.91 


dl 


0.29 +/- 


0.005 


0.08 +/- 


0.006 


0.32 +/- 


0.017 


0.081 +/- 


0.001 


0.046 +/- 


0.0011 


0.31 +/- 


0.007 


0.84 


dA 


0.31 +/- 


0.002 


0.11 +/- 


0.003 


0.32 +/- 


0.002 


0.066 +/- 


0.001 


0.049 +/- 


0.0009 


0.30 +/- 


0.004 


0.85 



a AII yields (Y) are [g/g]. Glucose (S), biomass (X), ethanol (E), glycerol (G), acetate (A), carbon dioxide (C02). 
b Strain abbreviations as in Figure 1. 



and easily folded, thereby minimally taxing the transla- 
tion and folding machinery. WA did increase two-fold 
in ATP consumption, most likely because a -amylase is 
10-fold larger and likely more difficult to fold and has 
more disulfide bond pairing possibilities. In the Ahacl 
background, folding efficiency is likely decreased due to 
ER dysfunction. With native secretion, dN did not 
require higher ATP maintenance consumption com- 
pared to WT. However, even the smaller, easier to fold 
IP resulted in ER stress that required significant ATP 
consumption compared to WT. dA, which was already 
stressed under WT, continued to show high ATP con- 
sumption. Despite the increased ATP consumption in dl 
and dA, little protein was secreted. 

Transcription factors controlling oxidative stress, amino 
acid salvaging, and expression repression are linked to 
secretory response 

Growth phase transcriptomics measurements were carried 
out to identify cellular processes that were activated under 
the stresses of HAC1 deletion and recombinant protein 
production. HAC1 deletion resulted in 339 significantly 
changed genes in the no recombinant protein case (WN 
vs. dN). HAC1 deletions in the insulin strain and a -amy- 
lase strain resulted in much larger cellular responses of 
1628 (WI vs. dl) and 1511 (WA vs. dA) significantly 
expressed genes, respectively. KAR2 (ER chaperone) 
expression was significantly reduced upon HAC1 deletion 
(| three-fold dN vs WTN, P=lx 10" 4 ) and the four yeast 
protein disulfide isomerases (PDI1, EUG1, MPD1, MPD2) 
reduced an average of 2.9-fold (P < 0.05). 

The effects of producing IP or a-amylase within a 
strain background (WT or HAC1) were not as pro- 
nounced as the effect of HAC1 deletion, 40 and 194 
genes were significantly changed in WI (compared to 
WN) and WA (compared to WN). Likwsise, 74 and 90 
genes were significantly changed for dl (compared to dN) 
and dA (compared to dN). 

To reduce the dimensionality of the data and identify 
putative TFs involved in protein secretion, the Reporter 
Transcription Factor algorithm was used [23]. TFs were 
scored by the modulation in expression level of genes 



that the TFs bind in the upstream region according to 
ChlP-chip data [24]. Therefore, the score is not indica- 
tive of change in the TF expression level itself, but of 
the genes under its influence. Reporter TF algorithm is 
useful, because although the statistical significance of an 
individual gene may not meet an arbitrary threshold, if 
several genes linked to the same TF have similar beha- 
vior, the likelihood of observing the group of genes is 
low, making TF identification very sensitive. Figure 3 
shows significant secretory process TFs shown to be 
involved in up- and down-regulating different cellular 
process under their control. Interestingly, different TFs 
were identified for the two different proteins. This is 
likely the combined effect of different protein size and 
number of disulfide bonds. A complete list of significant 
transcription factors is provided in Additional files 5 
and 6. 

In WT (Figure 3a), several TFs were activated by pro- 
tein secretion. Oxidative and osmotic stress pathway up- 
regulation was common to both proteins. Oxidative 
stress is likely caused by ROS that is formed when Erolp 
shuttles electrons to oxygen in disulfide bond formation 
[25]. Osmotic stress response, particular hypo-osmotic 
stress, strengthens the cell wall to counteract internal 
turgor pressure by changing the cell wall composition. 
This change in composition requires remodeling the 
secretory pathway by changing which components are 
trafficked to the cell wall [26] . Surprisingly, the Reporter 
TF algorithm found several Haclp-influenced genes 
down-regulated. Genes that Haclp binds from the ChlP- 
chip data that are significantly down-regulated are KEG1, 
MCD4, and ERJS. KEG1 and MCD4 genes are involved in 
glycan modifications and ERJS is a secondary ER chaper- 
one [27-29] . These genes may be influenced by other TFs 
not included in the ChlP-chip network. Genes known to 
be regulated by Haclp (KAR2 and EROl) were not signif- 
icantly changed upon secreting recombinant protein, 
indicating that there is not an actual Haclp response in 
the WT. 

Clear differences between large and small protein 
secretion emerge in WT. IP stimulated modification of 
the TGN through MCM1 and STE12. Overall expression 
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Figure 2 Secretory perturbations increase oxygen and ATP consumption, (a) Flux balance analysis of strains. TCA cycle is suppressed in 
high glucose. Bar graphs show flux in mmol/gDCW/h. Complete fluxes in Additional file 4. FBA calculations were performed using n = 3 
measurements with errors noted in Additional file 1 (Tables S1 and S2). The average coefficient of variance for these measurements is 11%. (b) 
Oxygen required for growth. The amount of oxygen consumed for each gram of cells (DCW) formed. Oxygen measurements are mean +/- s.e.m. 
(n = 3). (c) Specific maintenance ATP consumption as calculated by flux balance analysis (see Materials and Methods for details of calculation). 
FBA calculations were performed using n = 3 measurements with errors noted in Additional file 1 (Table S1 and S2). The average coefficient of 
variance for these measurements is 1 1%. Wild-type, no protein secretion (WN), wild-type insulin precursor secretion (Wl), wild-type a-amylase 
secretion (WA), Ahod no protein secretion (dN), Ahod insulin precursor secretion (dl), Ahacl a-amylase secretion (dA). Gram dry cell weight 
(gDCW). acetyl-CoA (Ac-CoA); fructose-6-phosphate (F6P); glyceraldehyde-3-phosphate (G3P); glucose-6-phosphate (G6P); pyruvate (PYR). 



is reduced by altering mRNA degradation pathways via 
STOl. a-amylase had a much larger effect on the cell, 
as compared to IP, as was implied by physiological para- 
meters of Figure 1 and number of altered genes. 



Additional oxidative and osmotic stress pathways were 
activated in WA, as well as a down-regulation in some 
amino acid synthesis pathways and overall reduction of 
transcription. 
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Figure 3 Transcription factors activated by recombinant protein secretion. Transcription factors were scored based on significantly 
changed genes in (a) wild-type strains, and (b) Ahacl. Venn diagram shows the number of secretory-related transcription factors activated in 
insulin precursor and a-amylase compared to no protein secretion. Table lists secretory-related transcription factors in small (insulin precursor) 
and large (a-amylase) protein secretion. Color coding indicates common secretory mechanisms as shown in diagram: modifying frans-Golgi 
network sorting (orange), oxidative stress (purple), amino acid metabolism (blue), and transcription and translation (green). 



In the Ahacl background (Figure 3b), many of the 
effects found in WA, have become common to both IP 
and a-amylase producing strains. HAC1 deletion clearly 
makes the cell more susceptible to recombinant secre- 
tion overload. Both insulin and a-amylase secretion 
cause considerable oxidative stress response and down- 
regulation of amino acid synthesis, including the general 
amino acid synthesis TF, Gcn4p. In dl, translational 
capacity repression is also employed (via Fhlp/Raplp) 
and adjustments in amino acid metabolism. dA shows a 
mix of up- and down-regulation of genes that are con- 
trolled by Haclp. Other TFs appear to be controlling 
these genes in the absence of HAC1. Some oxidative 
and osmotic stress pathways appear independent of 
HAC1. Skn7p and Cin5p were similarly activated in 
both WT and Ahacl. Oxidative and hypo-osmotic stress, 
while important for managing the secretory pathway, 
appears not to be directly managed through the UPR. 

Thermodynamic irreversibilities in redox reactions can 
explain increased oxidative stress in slow protein folding 
conditions 

The increases in oxidative stress, oxygen consumption, 
and reduced growth observed in the study can be 
explained by electron transfer in ER redox pathways. 
Disulfide bond formation has been established to con- 
sume oxygen and produce ROS (and thereby consume 
cellular resources to protect against the ROS) in stoi- 
chiometric quantities with the number of disulfide 
bonds formed [9]. When non-native disulfide linkages 
are formed, these linkages must be rearranged to the 



correct disulfide pairings for the native protein to be 
folded, a process called disulfide isomerization [30]. 

Disulfide isomerization involves (a) breaking the non- 
native bond by transferring electrons to the non-native 
bond creating a cysteine linkage with the PDI, and (b) 
creating a new disulfide linkage in the nascent protein 
by transferring the electrons to break the PDI-nascent 
protein linkage. By random pairing, the native disulfide 
bonds are found. 

Directionality in these redox reactions is determined 
by thermodynamic favorability through electron affinity 
of the potential disulfide bonds. Disulfide isomerization 
is redox neutral, not requiring electron donors or accep- 
tors. However, it does require each disulfide pairing to 
have a lower electron affinity than the next (non-native 
disulfide in folding protein < PDI-folding protein disul- 
fide < native disulfide in folding protein) to allow the 
electrons to transfer. Under slow folding conditions, PDI 
may hold the disulfide bond (oxidized state) for 
extended time because a native disulfide cannot be 
found, resulting in PDI being reduced by other moieties, 
likely GSH. 

Given the observations in our experiments, and the 
thermodynamic reasoning immediately above, we pro- 
pose a simple thermodynamic model of disulfide bond 
formation and breaking that explains increased oxidative 
stress, oxygen consumption, and reduced growth 
observed in our experiments. This model expands upon 
the mechanism by Cuozzo and Kaiser [13]. The thermo- 
dynamic model assumes there are PDI disulfide bonds 
that have electron affinities above and below the nascent 
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Figure 4 Proposed thermodynamic model predicts non-stoichiometric reactive oxygen species produced with incorrect disulfide bond 
formation, (a) In the model, forming and breaking an incorrect disulfide bond uses two protein disulfide isomerases (PDIs), one with electron 
affinity higher (PDI A ) and one lower (PDI B ) than the incorrect disulfide bond. In the formation phase, electrons are shuttled to molecular oxygen, 
resulting in ROS formation. In the breaking phase, electrons are passed from NADPH, through glutathione, to the protein. In both cases, 
electrons move along the electron affinity gradient. The net result is a futile cycle that is required to fix incorrect disulfide bonds, but expends 
redox energy, (b) The thermodynamic model predicts at fast folding rates near stoichiometric ROS is generated per disulfide bond formed. 
However, when folding rates are slow, the unfolded protein may go through many futile cycles, resulting in excess ROS. Glutathione (GSH), 
oxidized glutathione (GSSG), disulfide bond formation (DBF), disulfide bond breaking (DBB). 
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proteins disulfide bonds (Figure 4a). The disulfide is 
formed by the typical oxidation pathway (Figure 4a, 
green) catalyzed by high electron affinity PDI (called 
PDI A here). Instead of isomerization, the incorrect disul- 
fide is reduced by an electron donor with a low electron 
affinity (most likely a different PDI paralogue, called 
PDI B here) (Figure 4a, blue). The difference in electron 
affinity between the folding protein's cysteines and a 
specific PDFs cysteines can only allow the electrons to 
flow in one direction (toward the higher electron affinity 
cysteines) (Figure 4a). Therefore, a different PDI is 
required to form and break the incorrect disulfide bond. 
This futile cycle relies on a strong electron affinity gra- 
dient to complete an isomerization-like process. The net 
result of the futile cycle is GSH consumption and ROS 
production. This model implies that the ROS produced 
is not stoichiometrically linked to the number of disul- 
fide bonds formed, but varies by the number of futile 
cycles before the correct bond is formed. 

The metabolic and transcriptional data supports this 
model. Upon HAC1 deletion, ER chaperones (KAR2) and 
PDIs (PDI1, FUG1, MPD1, and MPD2) expression is 
reduced. This downregulation of ER chaperones and 
PDIs results in suppressed ER folding and disulfide bond 
formation in the Ahacl mutants. In the dN case, minimal 
oxidation stress is seen. However, when there is an 
increased demand for protein folding and disulfide bond 
formation, as is the case for dl and dA case, we see high 
oxygen consumption, ATP requirements, and many oxi- 
dative stress pathways being activated transcriptionally. 
Although both folding and disulfide bond formation is 
down, an imbalance toward faster disulfide bond forma- 
tion compared to folding will result in futile cycles. 
Therefore, this disulfide/folding imbalance acts as a cata- 
lyst for drastically increasing ROS production. 

Based on this thermodynamic model, the relative rates 
of protein folding and disulfide bond formation for nas- 
cent peptides have important consequences for oxidative 
stress (Figure 4b). When folding is faster than disulfide 
bond formation, ROS is produced in near one-to-one 
amounts with the disulfide bonds formed. Under these 
conditions, isomerization may be more efficient to resort 
incorrect disulfide bonds, as native structures with low 
electron affinity disulfide pairs are favored, and isomeri- 
zation does not produce ROS. However, when folding is 
slow compared to disulfide bond formation, as is the case 
when the protein folding machinery gets overloaded, the 
nascent peptides cycles through the futile redox cycle 
producing ROS and consuming GSH in excess to the 
final number of disulfide bonds formed. The physiologi- 
cal result of a high disulfide bond formation to ER folding 
rate is oxidative damage to a broad range of cellular pro- 
teins and consumption of reducing equivalents that 
could otherwise be used for anabolism. 



Discussion 

In this study, we have identified biological mechanisms 
related to protein synthesis and secretion by introducing 
perturbations to the cell, in the form of HAC1 deletion 
and different recombinant protein expression, and mea- 
suring the system level cellular responses, via transcrip- 
tomics and metabolic fluxes. These measurements, 
combined with data analysis algorithms, Reporter TF 
algorithm and FBA, were able to identify cellular adjust- 
ments in (a) overall expression level, (b) post-Golgi sort- 
ing, (c) amino acid biosynthesis and savaging, and (d) 
oxidative stress. These biological effects are a result of 
the combined influence of protein synthesis and traffick- 
ing through the secretory pathway. 

Overall transcription and translation were repressed in 
response to a -amylase expression (a larger protein) and 
in the Ahacl strains with any recombinant protein 
secretion. Repressing overall expression is a broad spec- 
trum response used to adjust the rates of all other cellu- 
lar processes to match the reduced folding capacity in 
the ER. Several mechanisms were used to alter overall 
expression: repressing mRNA synthesis, increasing 
mRNA degradation rates, and repressing protein transla- 
tion rates through reducing ribosome numbers. Specifi- 
cally, mRNA concentrations are lowered by decreasing 
RNA polymerase accessibility {HIR2), inhibiting tran- 
scriptional elongation (77/02), and controlling RNA 
degradation (STOl) [31,32]. Ribosome concentration, 
and thereby translation rates, can be reduced by the TFs 
Fhllp and Raplp which control expression of rRNA and 
ribosomal proteins [33]. This is seen in IP production in 
Ahacl strain, both by the reporter TFs (Figure 3) and 
by expression of ribosomal proteins (Additional file 7). 
In this context, extrachromosomal plasmids offer advan- 
tages over chromosomal expression. HIR2, whose 
mechanism is to silence the chromosome, would not 
affect extrachromosomal plasmids. Increased recombi- 
nant protein secretion would be accomplished by silen- 
cing native ER genes, while recombinant, plasmid-born 
gene would not be affected. 

Pronounced adjustments to the TGN were observed in 
the transcriptome in all conditions. TFs involved in phero- 
mone responses (STE12, MCM1, ASH1), invasive/pseudo- 
hyphal growth (STE12, MSN1, PHD1, RIM101), and 
osmotic stress (CIN5, SKN7, SKOl, YAP6, MSN1) were all 
identified by the Reporter TF algorithm and point to an 
underlying set of activities that are required to increase 
the traffic of secretory vesicles to the membrane. Invasive, 
pseudohyphal, and filamentous growth morphologies have 
a high surface to volume ratio and inherently require 
higher Golgi-to-cell membrane trafficking rates to supply 
cell membrane and cell wall components for growth. 
These altered morphologies can be activated through the 
filamentous and invasive response elements (FREs) [34] 
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bound by STE12 and used to regulate PHD1 [35]. HAC1 
deletion has been shown to cause filamentous growth [36]. 

Osmotic stress TFs are also responsible for affecting 
protein secretion, as the external cell wall must be 
strengthened in response to hypo-osmotic conditions, 
thereby requiring an efficient secretory pathway to ferry 
cell wall proteins [26]. MSN1 is known to induce starch 
degradation, requiring the actions necessary to secrete 
the appropriate enzymes through filamentous growth 
activation [37] . SKN7 has a dual role in invasive growth 
and osmotic stress [38]. Although osmotic stress TFs 
are commonly associated with the hyper-osmotic gly- 
cerol (HOG) pathway, Ypdlp can phosphorylate Skn7p, 
signaling the hypo-osmotic stress pathway [39]. Because 
there were no apparent hypo-osmotic conditions in this 
study, this indicates that these TFs are not directly con- 
trolled by osmotic conditions, but possibly through a 
secondary response to upregulation and increased secre- 
tion of cell wall proteins. 

TGN TFs and/or the genes they regulate are possible 
targets for increasing Golgi-to-cell membrane trafficking. 
In S. cerevisiae, recombinant protein intended for secre- 
tion has been found mis-trafficked to the vacuole. This 
has been shown for insulin and green fluorescent pro- 
tein secretion in yeast [40,41]. Proteins involved in vesi- 
cle trafficking, namely Slyl and Muncl8 have been 
found to increase recombinant secretory rates in Chi- 
nese hamster ovarian (CHO) and several mammalian 
cell lines [42,43]. It is likely that similar proteins are 
present in yeast and could be exploited for improving 
protein secretion. 

Significant alterations in amino acid metabolism were 
observed, particularly in the Ahacl strains. De novo 
amino acid synthesis {GCN4, BAS1, MET32, ARG81, 
RTG3) was suppressed. On the surface, this appears 
contradictory, as increased amino acid requirements 
should be observed with recombinant protein produc- 
tion. However, this decrease in de novo amino acid 
synthesis is accompanied by observed increases in 
scavenging mechanisms for amino acids (SNT2, CUP9, 
PUT3). High scavenging rates and decrease synthesis 
imply high protein degradation rates where the degraded 
proteins result in available amino acids for scavenging; 
reducing the need for newly synthesized amino acids. 
This is consistent with either ERAD, a process where 
proteins that are stalled in the ER are transported back 
into the cytoplasm for degradation by the proteosome, 
or vacuolar-localized protein degradation. In either case, 
the cell is expending energy on synthesizing proteins 
that are ultimately degraded. These effects appear in the 
strains that are the slowest growing with the highest 
ATP requirements (Figures lb and 2b). In these cases 
the ER folding capacity is likely saturated, resulting in 
ER holdup and ERAD. 



Oxidative stress TFs were also found in all conditions. 
Several were dual oxidative/osmotic stress TFs (CIN5, 
SKN7, SKOl), and others were dedicated to oxidative 
stress only (AFT2, YAP1). TFs were found in all three of 
the major oxidative stress signaling pathways, (a) the 
Hogl MAPK pathway (where SKOl is the DNA binding 
agent), (b) Slnl pathway (where SKN7 is the DNA bind- 
ing agent), and (c) YAP1 and CCS/5, which directly sense 
oxidative stress and bind DNA [44]. The cell's control 
machinery appears to have hard-wired oxidative stress 
responses to increased secretory demand, as oxidative/ 
hypo-osmotic pathways have a high degree of overlap, 
which is appropriate because increased secretion of cell 
wall proteins will result in higher oxidative stress. In 
particular, Skn7p, which has already been mentioned for 
its role in managing secretory pathway directly in an 
osmotic stress pathway, can also activate oxidative stress 
response genes [45]. 

Oxidative stress was pronounced with all secretory 
perturbations and has been identified in other studies to 
be associated with secretory stress [1,17]. Futile cycling 
may be the dominant disulfide resorting pathway when 
folding is limited. In previous studies, oxidative stress, 
induced by tunicamycin, a N-linked glycosylation inhibi- 
tor, increased with ER stress, despite no increase in the 
net disulfide bond formation demand [17]. The futile 
cycle does predict non-stoichiometric ROS formation, 
while isomerization does not. ROS can be formed at 
potentially limitless amounts through multiple rounds of 
disulfide formation and breaking. This will occur under 
conditions where the rate of folding is slow, a result of 
proteins that are specifically difficult to fold, or a result 
of the overall ER folding capacity being saturated. As 
well, futile cycling will increase as the number of avail- 
able cysteine residues available for disulfide bonding 
increase, as is the case for a-amylase, due to the 
extended amount of isomerization that may be needed 
to form the correct disulfide bonds. 

One implication of the proposed thermodynamic 
model is that PDI paralogues, or cysteines within a PDI, 
must exist at different electron affinities that are above 
and below the electron affinity of the protein to be 
folded. Although in vivo redox potentials of PDI cysteine 
pairs were not measured, from first principles it would 
appear highly likely that these PDIs would need different 
redox potentials to carry out isomerization. In Figure 4a, 
we assume that only PDIs interact with the folding pro- 
tein. This appears the case, as kinetic rates for direct 
glutathione oxidation/reduction are too slow to be phy- 
siologically relevant [9]. Electron affinity (and therefore 
redox potential) is broadly determined by the proximity 
of the two cysteines, with the proximity determined by 
the current structure of the protein [46]. Cysteines that 
are in the correct orientation will have a low electron 
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affinity and easily form disulfide bonds, while cysteines 
that are not in the correct orientation will have a high 
electron affinity and will have unstable disulfide bonds. 
Therefore, the electron affinity of a correctly folded/cor- 
rect disulfide bond would be lower than that of a mis- 
folded or incorrect disulfide bond. This difference in 
electron affinity may allow PDIs to selectively break dis- 
ulfides with high electron affinity (incorrect bonds), but 
not disulfide bonds with low affinity (correct bonds). 

The need for different PDIs to form or break disulfide 
bonds may explain the need for many PDI homologues 
in the ER, each with different structures, and therefore 
different electron affinities. These PDIs can only span a 
finite range of electron affinities, and there may be 
implications for proteins that have disulfide pairs with 
electron affinities higher than the highest PDI or lower 
than the lowest PDI. If no PDI has a lower electron affi- 
nity than an incorrect disulfide bond, then the disulfide 
bond cannot be broken and the protein is terminally 
misfolded. As well, a protein that has a native disulfide 
pairing with an electron affinity higher than any PDI 
cannot form a bond. This may be the case when recom- 
binant proteins are being processed in the ER. 

Futile cycling as a large potential ROS source has broad 
implications on the cell. Tu and Weissman predict Erolp- 
produced ROS that is one-to-one with disulfide bond for- 
mation could attribute approximately 25% of cellular ROS 
to the secretory pathway [1]. Therefore, even larger ROS 
production is likely if the futile cycle is the dominant dis- 
ulfide resorting pathway under folding stress. This also has 
implications on GSH and possibly NADPH availability, as 
it is doubly consumed (a) by the reduction of ROS and (b) 
directly in the futile cycle. The futile cycle limits reducing 
equivalents needed for anabolic processes, and may 
explain the reduced growth rates observed in folding 
stressed strains (WA, dl, and dA). 

In all, Figure 4b highlights that the relative rates of 
two processes, protein folding and disulfide bond forma- 
tion, must be kept in balance to avoid significant cellular 
stress. If disulfide bond formation is fast compared to 
folding, high futile cycle use will result in high ROS for- 
mation, NADPH loss, and high protein degradation as a 
result of ERAD. This scenario is observed in the Ahacl 
strains dl and dA. 

The engineering implications for protein secretion 
become much clearer with this understanding of protein 
folding to disulfide bond formation ratio. When overex- 
pressing a recombinant protein, an optimal expression 
must be found, where transcription is as high as possible 
without overloading the ER folding capacity and sending 
the cell into an oxidative stressed state. This optimal 
expression level will be different for different proteins, 
as protein folding rates will vary according to the pro- 
tein size and structure. We see this in comparing IP and 



a-amylase expression. The concept of an optimal 
expression has been identified heuristically, in the pre- 
sent study we identify the competing molecular effects 
that could define these phenomena [47]. This optimal 
expression ratio extends to recombinant proteins that 
do not have disulfide bonds. For recombinant proteins 
without disulfide bonds, recombinant protein folding in 
the ER will consume folding resources, thus slowing 
down folding rates. Although the recombinant protein 
has no disulfide bonds, many native proteins still require 
disulfide bonds. Because of this, the folding to disulfide 
bond formation ratio will be disturbed, resulting in simi- 
lar ROS stress. 

To maintain an optimal ratio, either protein folding 
rates must increase or oxidation rates decrease. Overex- 
pression of chaperones that increase folding capacity has 
successfully been used to increase protein secretion 
[6,48]. For particularly large or difficult to fold proteins 
this may not be adequate. A new approach may be to 
limit the oxidation rate of Erolp to slow down the first 
step of the futile cycle. This would be done in concert 
with repressing ERAD, as proteins would have long 
retention times in the ER. In this scenario, recombinant 
proteins would be slowly folded, albeit without high cel- 
lular stresses. This would result in longer overall process 
times, but may be required for difficult to fold proteins. 

Conclusion 

In this study, we identified post-Golgi vesicle sorting, 
high protein degradation rates, repressed overall expres- 
sion, and oxidative stress in response to +/- UPR strains 
secreting different sized recombinant protein. These 
processes were identified through scoring TFs and esti- 
mating alteration to the metabolic network. These 
observations imply our proposed futile cycling is the 
dominant disulfide resorting pathway in the ER and 
explains non-stoichiometric ROS formation seen in our 
study and elsewhere. The futile cycle model, producing 
ROS and consuming GSH, has a clear thermodynamic 
driving force compared to disulfide bond isomerization. 
If correct, futile cycling is likely the dominant mechan- 
ism under secretory stress. This interplay between pro- 
tein folding and futile cycling sheds light on a largely 
empirical understanding of engineering protein secretion 
and implies the relative rates of protein folding and dis- 
ulfide bond formation are critical to maintaining cellular 
homeostasis. This increased molecular understanding of 
the secretory pathway should allow for more insightful 
design of secretory engineering strategies. 

Methods 

Strains and media 

All experiments were performed in the background of 
CEN.PK 113-5D {MAT a SUC2 MAL2-8 C ura3-S2 1 P. 
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Kotter, Frankfurt, Germany) [49]. Genomic DNA from 
Y05650 (BY4741; Mat a; his3Dl; leu2D0; metl5D0; 
ura3D0; YFL031w::kanMX4, obtained from EURO- 
SCARF) was used as a template for the HAC1 knockout 
cassette. Standard molecular biology techniques were 
used [50] and all plasmids were maintained in Escheri- 
chia coli DH5a in Luria Bertani (LB) broth with 80 mg/L 
ampicillin. PCR primers are listed in Additional file 8. 

Cloning 

Genomic DNA was purified from Y05650 using Fast 
DNA Spin Kit for Soil (MP Biomedicals Solon, OH, 
USA). A 2.6 kb DNA fragment containing the genomic 
replacement of HAC1 with KanMX and flanking regions 
was amplified by PCR using primers KT007/KT008 
(Additional file 8). The 7MC7::kanMX4 fragment was 
integrated at the HAC1 loci of CEN.PK 113-5D by stan- 
dard yeast transformation [51] and selected on 200 mg/ 
L G418 to create the Ahacl strain. Correct integration 
was confirmed by PCR. 

DNA coding for an insulin precursor with a Yap3 pre- 
leader sequence and the TA57 pro-leader sequence and 
spacers as described [18] for correct secretory processing 
was synthesized with optimal codon usage for yeast and 
delivered on plasmid pUC57-Yap3Insulin (GenScript Co. 
Piscataway, NJ, USA) (Additional file 9 for sequence), a- 
Amylase DNA was amplified from Saccharomyces kluyveri 
YKM37 [52] using LZH018 and LZH039. The pre-pro-lea- 
der was amplified from pUC57-Yap3Insulin using primers 
LZH015 and LZH016. The pre-pro-leader was connected 
to the a-amylase by fusion PCR of the two segments 
together using primers LZH015 and LZH039 [53]. The 
pre-pro-insulin and pre-pro-amylase were cloned into the 
Spel/Sall or Spel/EcoRI sites of p426GPD, respectively, 
downstream of the constitutive GAPDH promoter [54], to 
create pYapIns and pYapAmy. Plasmids p426GPD, pYa- 
plns, and pYapAmy were transformed into CEN.PK 113- 
5D and Ahacl strains by standard methods [51]. 

Fermentor conditions 

Strains were grown in SD-2xSCAA [55], containing 20 g/L 
glucose, 6.7 g/L yeast nitrogen base minus amino acids 
(Formedium, Norfolk, UK), 2 g/L KH 2 P0 4 (pH = 6 by 
NaOH), 190 mg/L Arg, 108 mg/L Met, 52 mg/L Tyr, 290 
mg/L lie, 440 mg/L Lys, 200 mg/L Phe, 1260 mg/L Glu, 
400 mg/L Asp, 380 mg/L Val, 220 mg/L Thr, 130 mg/L 
Gly, 400 mg/L Leu, 40 mg/L Trp, 140 mg/L His, 1 g/L 
bovine serum albumin. Five hundred mL of medium was 
inoculated in a 1 L bioreactor (DasGip, Julich, Germany) 
at 30°C, 600 rpm agitation, 30 standard L/h air flow, pH 
controlled at 6 by KOH (2 M). Strains were inoculated to 
an A 60 o = 0.01 from late exponential phase cultures and 
Aeoo was measured throughout the cultivation. Dry cell 
weight (DCW) was measured by filtering 5 mL of culture 



broth through a 0.45 \im nitrocellulose filter and measur- 
ing the increased weight of the dry filter. Glucose, ethanol, 
glycerol, and acetate were measured using a Summit 
HPLC (Dionex, Thermo Scientific, Waltham, MA, USA) 
with an Aminex HPX-87H column (Bio-Rad, Hercules, 
CA, USA). Carbon dioxide and oxygen levels were mea- 
sured in the off-gas and dissolved oxygen was monitored. 
Transcriptome samples were taken after 5+ doublings at 
A 600 = 1.0-1.4. Triplicate fermentations were carried out 
for each strain. 

Protein quantification 

Insulin was measured by a modification of the assay by 
Snell et al. [56]. One mL of cell culture was centrifuged 
at 4000 x g for 4 min. Eight parts supernatant was 
added to one part 0.1 N HC1 and 5.5 (iM sodium azide 
and stored at 4°C until measurement. Insulin concentra- 
tion was determined by HPLC using a Luna 5 [i CI 8(2) 
(250 mm x 4.6 mm) (Phenomenex, Torrance, CA, USA) 
column and gradient-based elution. Buffer A contained 
68 mM phosphoric acid, 0.2 M sodiumsulphate and 10% 
(w/v) acetonitrile in water, and Buffer B contained 50% 
acetonitrile in water. HPLC was run with 25 [iL injec- 
tions at 1 mL/min and 50°C. Gradient protocol: 20% B 
for 10 min. Linear gradient from 20% B to 60% B over 
10 min. Hold at 60% B for 5 min and then to 20% B for 
3 min to re-equilibrate for next sample. Insulin stan- 
dards eluted at 22.6 min and insulin precursor at 20.0 
min. HPLC peaks were verified to be the correct protein 
by SDS-PAGE. Human insulin was used as a standard 
(Sigma, St. Louis, MO, USA). 

a-amylase concentration was calculated from enzyme 
activity, a-amylase activity was measured using the Cer- 
alpha kit (Megazyme K-CERA, Wicklow, Ireland) using 
a-amylase from Aspergillus oryzae (Sigma, St. Louis, 
MO, USA) as a standard. This conversion was calculated 
using a 1.79 U/mg (weight includes salts and purified 
protein) standard from Sigma using the Protein 80 chip 
on the Bioanalyzer (Agilent, Santa Clara, CA, USA). By 
this, a-amylase was found to be 0.0257 g a-amylase/g 
total. a-amylase activity was then converted to mass 
using 70 U/mg a-amylase protein. 

Transcriptome analysis 

Samples for microarray were taken as described pre- 
viously and stored at -80°C until processing [57]. RNA 
was isolated using the RNeasy Minikit (Qiagen, Valen- 
cia, CA, USA). Cells were lysed in RNeasy RLT buffer 
using Lysing Matrix C (MP Biomedicals Solon, OH, 
USA) in a Fast Prep 24 (MP Biomedicals Solon, OH, 
USA) as follows: 20 s at speed 6, 1 min at 4°C, 20 s at 
speed 6. RNA was processed to aRNA using the Gene- 
chip 3' IVT Express Kit (Affymetrix, Santa Clara, CA, 
USA) and hybridized/scanned on the Yeast Genome 2.0 
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Array (Affymetrix, Santa Clara, CA, USA) following 
commercial protocols to create CEL files. 

Images were analyzed using R 2.10.1 statistical soft- 
ware and the 'affy' and 'limma' packages as described 
previously [58]. Briefly, background normalization was 
carried out using robust multi-array (RMA) average 
method with perfect match (PM) probes only. Interchip 
normalization used the qspline algorithm with median 
polish summary method. Statistical analysis was carried 
out by comparison of triplicate bioreactor measurements 
for each strain. Emperical Bayesian statistics were used 
to moderate standard errors within each gene and Ben- 
jamini-Hochberg's method to adjust for multiple testing. 
Microarray data was submitted to the GEO database 
and have accession number GSE27062 (see http://www. 
ncbi.nlm.nih.gov/geo/query/acc.cgi? 
token=dpyzfywysoqecbk&acc=GSE27062). 

Reporter transcription factor analysis 

Transcription factor activity was scored using the 
Reporter Effector algorithm [23]. Transcription factor- 
DNA interactions were gathered from ChlP-chip with P 
< 0.001 [24]. Significant interactions were found for 176 
transcription factors regulating 3,796 genes for a total of 
10,849 unique interactions. Gene P values from compar- 
ing different strains were used to score transcription fac- 
tors that were known to bind to the upstream DNA. 
Transcription factors with P < 0.05 of being activated 
between conditions are reported. 

Flux balance analysis 

Estimates of intracellular reaction rates were performed 
using measured exchange fluxes of glucose, ethanol, 
acetate, glycerol, and carbon dioxide. Model-based error 
correction was used to close carbon and electron bal- 
ances [59]. Flux balance analysis was carried out using a 
85 reaction model of yeast central carbon metabolism 
and biomass yield were used [21]. Additional file 4 con- 
tains the complete results of the analysis which are used 
to estimate ATP consumption in the different strains. 

Additional material 



Additional file 1: Measured exchange fluxes in strains. Measured 
metabolite exchange fluxes for strains used in this study. 

Additional file 2: Final glycerol concentration of WT and Abaci 
strains. Measured glycerol titers at end of fermentation for strains used 
in this study. 

Additional file 3: Estimated exchange fluxes. Metabolite exchange 
fluxes as estimated by error-correction algorithm for strains in this study. 

Additional file 4: Intracellular fluxes for metabolic network. Flux 
balance analysis estimates of internal fluxes for strains in thus study. 

Additional file 5: Reporter TFs for WT protein secretion. Transcription 
factors activated by recombinant protein secretion in wild-type 
background. 



Additional file 6: Reporter TFs for Ahacl protein secretion. 

Transcription factors activated by recombinant protein secretion in Ahacl 
background. 

Additional file 7: Expression profiles for ribosomal proteins. mRNA 
concentrations for yeast ribosomal proteins as determined by DNA 
microarray. 

Additional file 8: Oligonucleotides used in this study. PCR primers 
used for cloning and validation. 

Additional file 9: Synthesized insulin precursor DNA sequence. DNA 

sequence for insulin precursor used in this study. 
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